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Abstract 

We present a measurement of the meson lifetime in the decay mode J /ipe'^Ue using 

the CDF II detector at the Fermilab Tevatron Cohider. From a sample of about 360 pb~^ of pp 
collisions at ^/s = 1.96 TeV, we reconstruct J/ipe'^ pairs with invariant mass in the kinematically 
allowed range 4 < Mj/^g < 6GeV/c^. A fit to the decay- length distribution of 238 signal events 

yields a measured B^ meson lifetime of 0.463 ^q'qq^ (stat) ± 0.036 (syst) ps. 

PACS numbers: 13.20.He, 14.40.Nd 



The meson is the only known meson consisting of two heavy quarks of different 
flavor: a charm quark and a bottom antiquark. It provides a unique test of heavy-quark 
dynamics, since the bound state can be treated using the same non-relativistic expansion 
that successfully describes both cc and bb families. However, unlike cc and bb states, the 
meson decays only via weak interactions, thus having a measurable lifetime. The lifetime 
of the B^ meson is expected to be about two to three times smaller than the B~^ meson 
lifetime if one assumes three major decay subprocesses b quark decay with c quark as 

spectator, c quark decay with b quark as spectator, and be annihilation decays. In the B~^ 
meson case, the dominant decay subprocess is the b quark decay with u quark as spectator. 
An early measurement from CDF found a B^ lifetime consistent with predictions. More 
precise measurements will determine the relative importance of the three decay subprocesses 
and provide insight into the strong dynamics of heavy quarks. Here, we report a new B^ 
meson lifetime measurement using the decay mode B'^ J/ipe^Uf., where charge-conjugate 
modes are implied, from a sample of 360 pb^^ collected during 2002-2004 with the CDF II 
detector |^ at the Fermilab Tevatron Collider at a center of mass energy of 1.96 TeV. The 
new measurement is about 2.5 times more precise than the previous one. 

The B^ J/ipe'^Ve reconstruction starts with J/ip ^ /^"'"yU" candidates selected based 
upon a two-muon topology by the CDF trigger system j^. The J/ifj candidates are further 
purified during offline reconstruction by vertex constraining the fi^fi' pairs and by selecting 
the pairs with with momentum transverse to the beam line pTj/^>'i GeV/c. Then, each J 
candidate with a reconstructed mass within 50 MeV/c^ of its nominal value is combined with 
an electron to form a B^ candidate. 

Electron identification uses both specific ionization {dE / dx) information from the Central 
Outer Tracker (COT) and calorimeter shower information from the Central ElectroMagnetic 
calorimeter (CEM). The logarithm of the ratio of the measured dE / dx value from a charged 
particle to that expected for an electron, Zf, = \n{dE / dx) — In (c^E/fix) predict, is compared 
to its standard deviation az^ ■ The expected dE / dx and cr^^ are functions of the particle 
charge, momentum, and the multiplicity of associated COT hits. Electron candidates are 
required to have Z^/az^ > —1.3 to reject hadrons (tc/K/p) while remaining efficient for 
true electrons. Samples of electrons, pions, kaons, and protons selected from collision data 
are used to determine the dE / dx identiflcation efficiencies listed in Table HI These control 
samples come from photon conversions, 7 e~^e~, and from hadron decays, K'^ —* tt~^tt~, 



7r+, and hP —>■ pn^ . The hadrons surviving the Z(,/crz^ selection are mainly 
pions, which are rejected using calorimeter shower shape information. 

The calorimeter shower shape of a charged particle with prp>2GeV/c is obtained by 
extrapolating its track reconstructed in the COT into the calorimeter to match shower 
clusters there The probabilities for a particle to have a shower shape consistent with 
being an electron or a hadron are calculated using the distributions of shower energy and 
shower cluster profiles for the electron and hadron samples described above. We first define 
the probability for a charged particle to be an electron based on shower shape by the ratio 
between its probability to be an electron and the sum of its probabilities to be an electron 
or a hadron. We then obtain a cumulative probability distribution of the ratio using the 
electron sample and impose a selection on electron candidates at a 70% probability value. To 
calculate the average probability for hadrons to pass this requirement, as listed in Table HI 
the control samples of 7c/K/p particles are mixed using fractions predicted by a PYTHIA 
Monte Carlo (MC) simulation of B —>■ J/ipX events [6]. In addition, electrons found to 
originate from photon conversion, 7 e~^e~, are removed from consideration as J/ipe^ 
candidates P]. Overall, the electron identification using combined clE / dx and calorimeter 
information has an efficiency of 60% for true electrons, while hadrons have a probability to 
pass the selection lower than 0.15%. 

TABLE I: Electron identification efficiencies (%) as functions of particle pT (GeV/c) using dE/dx 
and calorimeter (Cal) for electron and hadrons. The dE/dx results are averages of positively and 
negatively charged particles. The calorimeter results for hadrons {h) are the weighted averages of 
vr, K and p. The CO livers ion- fin ding efficiency Cconv 

is also listed. 



PT 


2-3 


3-4 


4-5 


5-6 


> 6 


dE/dx:e 


91.3 ±0.1 


91.4 ±0.2 


90.7 ±0.2 


90.5 ±0.3 


89.5 ±0.2 


dE/dx:n 


16.4 ±0.1 


23.8 ±0.1 


32.0 ±0.1 


39.0 ±0.2 


49.4 ±0.2 


dE/dx:K 


2.09 ± 0.09 


2.35 ±0.07 


3.29 ±0.09 


4.4 ±0.1 


9.5 ±0.2 


dE/dx:p 


3.38 ± 0.03 


2.14 ±0.04 


2.54 ±0.07 


3.0 ±0.1 


4.9 ±0.2 


Cal:e+ 


68.5 ±0.3 


68.8 ±0.5 


68.9 ± 0.8 


68.9 ± 1.1 


67.6 ± 1.0 


Cal:e~ 


67.9 ±0.4 


69.5 ±0.5 


69.6 ± 0.7 


68.1 ± 1.2 


68.5 ±0.9 


Cai.h+ 


0.77 ±0.04 


0.37 ±0.04 


0.37 ±0.03 


0.29 ± 0.04 


0.13 ±0.04 


Cal.h- 


0.64 ± 0.04 


0.37 ±0.02 


0.27 ±0.03 


0.25 ± 0.04 


0.21 ±0.04 


^conv 


49.8 ± 1.4 


55.0 ±2.2 


56.5 ± 3.3 


61.5 ±4.5 


69.2 ±3.4 



A candidate is a J/ipe~^ pair with transverse momentum ptj,^^>5 GeV/c and invariant 



mass 4<Mj/^e<6 GeV/c^. The upper bound on Mjj^f. is simply the kinematic hmit from 
the mass of meson 1\. The lower bound is set higher than the kinematic limit of 
Mjj^ to reduce the background from 5+ semileptonic decays other than the exclusive decay 
5+ J/ipe'^i'e to a few percent . The opening angle between the J/ip and electron 
momenta in the transverse plane must be within 90° to reduce generic hh background that 
produces a, J/ip and an electron from different h hadrons. Finally, the tracks of the three 
daughter particles, /x"*", and e"*", are fit to a common vertex and the 5+ decay length 
in the transverse plane, Lxy, is calculated as the projection of the displacement of the B'^ 
vertex from the primary vertex onto the momentum of the J/ipe^ system. The primary 
vertex position is obtained from run-by-run averages using samples of prompt tracks. 

Before making a lifetime measurement, we first establish the 5^ signal in the J/ipe^ 
pairs. The background pairs from prompt decays are removed by imposing a selection 
of LxylcrL^y> 3. The Mjj^^ distribution of J/ipe^ pairs with Lxy/aL^y> 3 is shown in 
FigHJ Within the 4<Mj/^e<6 GeV/c^ window, 203 candidates are found with a total of 
background 88 ± 14 as listed in Table |nj The background is classified into four groups: 
(i) background with a false J/V", background with a correctly identified J /if) but a wrong 
electron candidate which can either be (ii)a mis-identified hadron (false electron) or (iii) 
an electron from photon to e^e~ conversion, and (iv) background from decays of other 
B hadrons (66). The number of false J/ip background is estimated using pairs with 

invariant mass outside the 50 MeV/ window. To estimate the false electron background, we 
use a sample of J / ■j/'-track pairs passing the same selection as the J/ ipe^ pairs, including the 
dE I dx requirement and the requirement to point to the GEM fiducial region, but without any 
selection based on calorimeter information. The size of the contribution is estimated from a 
weighted counting of the J/ tp-tiack pairs with the weights taken as the averaged probabilities 
in Table H] for hadrons to pass the electron selection using the calorimeter. We derive 
the residual conversion electron contribution from the rate of identified photon conversions 
together with the conversion-finding efficiency listed in Table H] The conversion-finding 
efficiency, defined as the fraction of the identified electrons with their conversion partners in 
the kinematic acceptance of the GDF detector, is estimated from a MG simulation. Finally, 
the contribution from decays resulting from 66 production is estimated using a PYTHIA 
MG sample with relative rates of fiavor creation, fiavor excitation, and gluon splitting tuned 
to the Tevatron data 0, 0]. The number of B^ signal events is found to be 115 ± 16 (stat) ± 



14 (syst). For comparison, there are 2872 ± 59 J/tpK^ events in the data sample 

corresponding to the same integrated luminosity. The selection of B'^ is the same as for B^ 
except for electron identification on the K^. We find the production rate of B'^ relative to 
that of B+, [(r{B+) ■ B{B+ J /%l)Xe+Ve)]l\a{B+) ■ B{B+ J/ipK+)], in the kinematic 
range pr > 4GeV/c and rapidity \y\ < 1 7] to be 0.282 ± 0.038 {stat) ± 0.035 {syst) ± 
0.065 {acceptance) . The first error is statistical; the second covers the systematic uncertainty 
of B^ signal excess counting, and the third pertains to the estimated detector acceptance 
ratio correction Ab+ = 4.42 ± 1.02 from MC simulation where the B^ px spectrum, its 
lifetime values and decay modes are the major sources of uncertainty. This new production 
ratio result agrees with the earlier CDF measurement j3|. 
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FIG. 1: Mj/^e distribution of J/V'e"'" pairs (dots with error bars) with Lxy/crL^y>3 together with 
expected shape (dashed line) from a sum of a MC signal and estimated background (solid lines) 
which is mostly from 6-hadron decays. The false J/ip is already removed in the distribution. 



TABLE II: Numbers of J/tpe'^ pairs and estimated backgrounds. The error listed is the sum of 
statistical and systematic errors. 





All L^y 


Lxy/UL^^ > 3 


false J/-)/) 


164.0 ±9.1 


24.5 ±3.5 


false electron 


110.2 ± 19.0 


15.4 ± 2.5 


conversion electron 


67.4 ±34.8 


14.5 ± 7.8 


bb 


63.0 ±18.5 


33.6 ± 11.4 


prompt decay 


141.7 ±32.0 




total background 


545 ± 55 


88 ± 14 


observed J/tpe~^ pairs 


783 


203 



Having established a clear 5+ signal in J / ipe'^ combinations, we measure the B'^ meson 
lifetime in a larger sample of 783 events, selected with the same criteria as above, but 
without the L^y selection. We estimate the net signal excess in this sample to be 238. The 
background sources and their contributions are listed in Table IH] Contributions from false 
electrons, conversion electrons, 66, and false J/ip are estimated as described earlier. The 
number of additional prompt-decay events is extracted directly from the lifetime fit. 

The Lorentz-invariant proper decay time of a B'^ event is its decay length L^y with 
a Lorentz boost /?7 = ptb^/^b^ in the transverse plane, where Mb^ is the mass of B'^ 
meson and p^c. its transverse momentum. The -B^ mass is assumed to have a value of 
6.271 GeV/c^ Because of the missing neutrino, we cannot directly calculate the boost 
factor from the lab system to the 5^ rest frame. We can however calculate the B'^ lifetime 
in the center of mass frame of the J/ipe^ pair, ct' = LxyMj/^^/PTj/^^, which provides the 
best estimator of the B^ proper decay time in the absence of the neutrino momentum. The 
true proper i?+ decay time is given by ct = LxyPTB^I = K ■ ct', where K is the residual 
correction factor depending on the momenta of the missing neutrino. In MC simulations, 
K can be calculated K = {ptj/^^/ptbJ ' {^Sal ^J/^e) ■ cos a, where a is the angle between 
the vectors of ptj/^^ and Ptb^- The K distributions, as shown in FigEl always include K=l, 
and decrease in width as Mjj^f. gets closer to the value of B"^ mass when the neutrino's 
momentum is minimal. 

An unbinned maximum-likelihood fit 0, 0] is used to extract the B^ meson lifetime. 
In the fit, t! and its event-by-event error, cr^/ = ai^yMji^e/ {cpTj/^J^ are the input vari- 
ables. The likelihood function has the form [9| T{t',at') = (1 — J2i fbj^sit' , (Jti)Vs{(Jt') + 
J2i hi^biit' , ati)Vh.X^t') , where J-'sit' , af) is the lifetime Probability Density Function (PDF) 
for pure B'^ signal, /;,. and J-'bXt' , erf) are fractions of the five background contributions 
and their lifetime PDFs, and Vs{<yt') and VbX'^t') are the PDFs of the at' for signal and 
backgrounds. The lifetime PDF for the B^ signal is an exponential lifetime distribution 
convoluted with the K distribution and a Gaussian resolution function. The prompt back- 
ground lifetime PDF is assumed to have zero lifetime with a Gaussian resolution function. 
The PDFs for other backgrounds are described by a sum of a Gaussian distribution centered 
at zero and two pairs of positive and negative exponential lifetime functions with no K 
correction applied. The initial parameters for these background PDFs are obtained from fits 
to the background samples as shown in Fig. 121 The obtained results are used to constrain 



the corresponding parameters in the final hfetime fit. The constraints are imposed by mul- 
tiplying the likelihood function with Gaussian functions of the appropriate mean and width. 
Similarly, the background fractions fb. are also constrained in the B'^ meson lifetime fit to 
the estimated values in Table ITTl 



I" 



1- 




(a) 




5- f 


1 — ^ — I 



















0.05 



0.5 



1.5 




K distribution 



FIG. 2: Distribution of K for Mj/^^ within (a) 4-4.5, (b) 4.5-5.0, (c) 5.0-5.5 and (d) 5.5-6.0 GeV/c^. 
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FIG. 3: ct' PDF obtained from fits to background samples (points with error bars) for (a) false 
J/ip (solid lines), (b) bb (solid lines), (c) false electron (hatched area) and (d) conversion electron 
(hatched area). In (c) and (d), the solid lines are the sum of false J ftp (shaded area) and false or 
conversion electron. The false J ftp fraction and lifetime shape in (c) and (d) are constrained to 
that obtained from sideband events. 



In FiglH the ct' distribution from the 783 candidates is shown with the fit result 
superimposed. We find ctb^ = 139^20 A*^- The sources of systematic uncertainty on the 
lifetime fit are now considered, and their magnitudes are estimated. The effect of the K 
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FIG. 4: ct' distribution from 783 B'^ candidates. The points with error bars are data points and the 
sohd hne is the fit result. The dashed Une is the B'^ signal and the dotted line is the background. 

distribution uncertainty on the B'^ lifetime fit is estimated using alternative K distributions 
obtained from MC simulations with different pt spectrum, mass and lifetime values, and 
-B^ meson decay modes. The px spectrum is changed from that derived using a theoretical 
calculation 3] to that of the inclusive decay B ^ J / ipX [4] . The mass and lifetime in the 
MC calculation are varied in the ranges Mb^ = 6.2-6.4GeV/c^ and r^^ = 0.4-0.7ps 
The B^ decay considered in the MC simulation is varied from the exclusive B^ J /ipe^Ve 
alone to that of an inclusive decay table predicted in Ref. We found the change on 
the B'^ lifetime fit as Acr^^ = ±2.8 /im. The uncertainty related to background lifetime 
shapes is estimated from investigating the pt dependence of the electron identification and 
the conversion- finding efficiencies, from using false J/ip events from different samples with 
or without an electron nearby, and from changing fractions of hb events originating from the 
three main production mechanisms according to a study using CDF data js| . The estimated 
effect is ±9.2 /im. The uncertainty related to the L^y calculation and its error distribution is 
found to be ±4.7 /xm from the uncertainty in the silicon detector alignment, by using an al- 
ternative functional form describing the L^y resolution that includes an additional Gaussian 
and symmetric exponential tails, and by using alternative decay-length error distributions. 
The fitting procedure is also checked using MC ct' distributions similar to that of the 5+ can- 
didates and there is no bias found. Adding all the estimated systematic errors in quadrature 
we find ctb, = 139 tfo (stat) ± 11 (syst) fim or tb, = 0.463 t^oml {stat) ± 0.036 {syst) ps. 

In conclusion, from an unbinned maximum-likelihood fit to the decay-length distri- 
bution of 238 signal events of B'^ — > J/ipe^u^, the B'^ meson lifetime is found to be 



0.463 

-o'ml (stat) ± 0.036 (syst) ps which is about one-third of the 5+ meson hfetime 
of 1.671 ± 0.018 ps [7J. This agrees with theoretical models in which all the three 

major decay subprocesses, the two spectator processes (6-quark and c-quark) and the be 
annihilation, play important roles in the decays. 
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